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ABSTRACT 

Context. We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence 
on the illuminating UV field. 

Aims. Using Herschel-KWl observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the 
prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. 

Methods. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinc- 
tion by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically 
the intensity and line profile of the 158 /jm (1901 GHz) [Cn] line measured by HIFI and provide information on the emitting gas. 
Results. We show that both the [C n] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions lo- 
cated in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement 
between the calculated and observed [C n] intensities. 

Conclusions. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include 
analysis of the main coolant line [O i] and use of a new PDR model that includes PAH-related species. 
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1. Introduction 

One main goal of the Guarantee Time Key Programme "Warm 
and dense interstell ar medium" (WADI) o f the HIFI hetero- 
dyne spectrometer dde Graauw et al ] I2010h onboard Herschel 
dPilbratt et al.ll2010T) is to investigate the physics and chemistry 
of the gas and dust in dense photon-dominated regions (PDRs), 
as well as their dependence on the illuminating UV field. As part 
of this programme, we observed a prototype PDR, NGC 7023. 
The region is illuminated by the B2Ve HD 200775 [RA(2000) 
= 21h01m36.9s ; Dec(2000) = +68°09'47.8"], and has been 
shaped by the star formation process leading to the formation 
of a cavity. NGC 7023 has been widely studied at many wave- 
lengths. It has been shown that this region hosts structures at dif- 
ferent gas densities: n# ~ 100 cm -3 in the cavity, ~ 10 4 cm~ 3 in 
the PDRs that are located north-west (NW), south (S) and east, 
and 10 5 - 10 6 cm ~ 3 in dense filaments and clumps that are ob- 
served in the mm dFuente et alj 1996t Gerin et al [1991 and ref- 
erences therein) and near-IR ( Lemaire et al. 1996 ; Martini et alj 
119971) . 

NGC 7023 has been mapped by the instruments PACS and 
SPIRE of Herschel to study the emission of large cold grains 
dAbergel et alj|2010l) . We present here some observations of the 
gas at the surface of this nebula, taking advantage of the very 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia with impor- 
tant participation from NASA. 



high spectral resolution of HIFI. By combining these observa- 
tions with previous mid-IR observations, we study the geometry 
and energetics of the NW and S PDRs. 



2. Observations and results 

2.1. HIFI observations 

The HIFI observations presented here consist in (offset positions 
are relative to the star; see Fig.Q]): 

- single pointing using frequency switch mode towards the NW 
PDR in bands la and 3b. In band la, the frequency ranges 
covered by the Wide Band Spectrometer (WBS) were [535- 
539] GHz (LSB) and [547-551] GHz (USB). In band 3b, covered 
ranges were [879-883] GHz (LSB) and [891-895] GHz (USB). 
The observed position (Aa = -25", AS = +38", called H2 peak) 
corresponds to the peak intensity of the H2 ro-vibrational em is- 
sion associated to the near-IR filaments dLemaire et alJI 1 996h . 
-On-the-fly (OTF) mapping of the [C 11] 1901 GHz emission line 
in band 7b. The line was covered by both the WBS and the 
high-resolution spectrometer (HRS) in USB. Two cuts were per- 
formed: a cut from the star to the NW PDR (star-NW cut) and 
a south-north cut (hereafter S-NW cut) covering the NW PDR, 
the cavity and the S PDR (Aa = -47", -85 < AS < +60", see 
Fig.[TJ. The pixel size after regriding is 6.5". 

All these observations include an OFF reference position in 
the western lobe of the cavity (Aa = -144", AS = -47" ). 
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Table 1. Summary of the HIFI data 



Transition 


v line 


Position 


Beam 




FWHM T 


T , 


Area 




[GHz] 




size ["] 


[km s~'] 


[km s~'] 


[K] 


[ Kkms~'] 


HCO + (6-5) 


535.062 


H 2 peak 


41 


2.3 ±0.3 


1.4 ±0.3 


0.22 


0.33 ± 0.07 


13 CO (5-4) 


550.926 


H2 peak 


41 


2.2 ±0.3 


1.2 ±0.3 


5.12 


6.5 ± 1.5 


13 CO (8-7) 


881.272 


H 2 peak 


24 


2.2 ±0.2 


0.9 ±0.2 


4.15 


4.0 ±0.8 


C 18 (5-4) 


548.831 


H 2 peak 


41 


2.1 ±0.3 


0.9 ±0.3 


1.12 


1.1 ±0.4 


[Cn] 


1900.537 


Ht peak 


11 


2.7 ±0.1 


3.4 ±0.2 


23.6 


85.6 ±5.0 


[Cn] 


1900.537 


S-NW cut / North 


11 


2.5 ±0.1 


2.4 ±0.2 


17.9 


45.8 ±3.8 


[Cn] 


1900.537 


S-NW cut / Cavity 


11 


4.0 ±0.1 


2.2 ±0.2 


6.89 


16.2 ± 1.4 


[Cn] 


1900.537 


S-NW cut / South 


11 


2.5 ±0.1 


2.4 ±0.2 


13.3 


34.0 ± 2.8 



For [C 11] , we report the FWHM of the Gaussian profile of equivalent area and peak intensity. 
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Fig. 1. The NW and S PDRs of NGC 7023 observed by Spitzer- 
IRAC at 8 pm (red) and 3.6 pm (green). The white circle repre- 
sents the HIFI beam at 535 GHz (41") towards the H 2 peak. The 
dotted lines show the cuts that are studied in the [C 11] emission 
line at 158//m with a beam of 11", whereas white crosses in- 
dicate the specific positions reported in Fig. [2] and Table [TJ The 
star position is shown with a black cross. 

Data was reduced with HIPE 3.0 (IOttll2010h on level-2 data pro- 
duced with the standard pipeline. For the pointed observations in 
bands la and 3b, manual steps consisted in stitching sub-bands, 
baseline removal, and correction for main beam efficiency (77,,,/, 
= 0.71). The [C 11] WBS spectra required defringing, which was 
performed with standard HIPE tasks, and the best data quality 
was produced with the subtraction of two sinusoidal fringes. To 
verify the biases introduced by the fringe removal, we compared 
the WBS and HRS profiles, which showed good agreement both 
in profile and absolute intensity except for the weakest lines 
(T* <4K). 

2.2. Gas kinematics 

Figure|2](left panel) shows the 13 CO 5-4, 13 CO 8-7, C ls O 5-4, 
and HCO + 6-5 lines observed by HIFI towards the H2 peak. All 
the lines have a central velocity of about 2.2 km s , comparable 
to previous ground- based observations in several molecular lines 
dFuente et alJll993l) . Figure [2] (right panel) shows the [Cn] line 
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Fig. 2. Left: 13 CO, C ls O, and HCO + high-J transitions observed 
with HIFI toward the H2 peak. Right: Examples of [C n] emis- 
sion profiles at the H2 peak and at the three positions shown 
by crosses on Fig.Q] north (Aa = -47", A6 = +45"), cavity 
(Aa = -47", A5=-20"), and south (Aa = -47", A6 = -80"), 
with offsets relative to the star position. In both figures, the ver- 
tical dotted line indicates v/ sr = 2.2 km s _1 . 



profiles at the H2 peak and at different positions along the S-NW 
cut. The line is much broader than molecular lines and its profile 
shows a complex multi-component structure. The observations 
towards the PDRs show that the emission peaks at intermediate 
velocities (v/„. ~1. 8-2. 8km s ), which have already been ob- 
served toward s the NW PDR in several molecular tracers (see, 
for example, iFuente et al.lll996T) . There is also a contribution 
from higher velocity components (v/„ ~4km s _1 ), which domi- 
nate the emission in the cavity. 

2.3. Emission from very small dust particles and C + 

The 158 yum (1901 GHz) [Cn] and 63 //m [ O i] lines are the major 
coola nts of the gas at the surface of PDRs dHollenbach & Tielensl 
[19991) . In these regions, photoelectric effect dominates the heat- 
ing, while H2 formation provides a minor contribution. Since 
the smallest dust particles, polycyclic aromatic hydrocarbons 
(PAHs) and very sm all grains (VSGs) contribute to a large frac - 
tion to this process (iBakes & Tielenslll994tlHabart et aljbool . 
and these particles emit in the mid-IR most of the energy they 
absorb in the UV, we expect that the mid-IR and [C n] emissions 
arise in the same regions. 
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Fig. 3. Comparison between the [Cn] 158jum line flux (solid 
line) measured with HIFI at a beam size of 11" and the aro- 
matic IR band (AIB) flux (5.5-14/im) along the star-NW (a) and 
S-NW (b) cuts. The error bars for [C n] are computed at one- 
sigma level. The AIB flux is determined with a fit of the mid-IR 
spectra using the three PAH-related populations shown in Fig|5j 
filled diamonds are Spitzer data (1.8"pix _1 ), and open diamonds 
are ISOCAM data ^"pix" 1 ). 
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Fig. 4. Comparison between the [Cn] 158 //m line flux (solid 
line) measured with HIFI at a beam size of 1 1" and the column 
density N(H) along the star-NW (a) and S-NW (b) cuts. N(H) 
was derived from both the AIB flux (diamonds) and the mid-IR 
dust extinction (open circles); filled diamonds and open circles 
are Spitzer data (1.8"pix~'), and open diamonds are ISOCAM 
data (6"pix~'). 



We used mid-IR spectro-imagery data of the NW PDR of 
NGC 7023 that were obtained in the 5.5 -14.5 range wit h 
the Infrared Spectrograph onboard Spitzer dWerner et alj|2004l) . 
For t he S PDR, we used IS OCAM highly-processed data prod- 
ucts dBoulanger et al.|[l996l) from the ISO data archive. To anal- 
yse the mid-IR spect ra, we followed the m ethod explained in 
iRapacioli etal] (120051) and lBerne etall (120071) . in which the mid- 
IR emission is decomposed into three aromatic IR band (AIB) 
spectra whose carriers are neutral PAHs (PAH ), cationic PAHs 
(PAH + ), and evaporating VSGs. The fitting procedure was re- 
cently improved by inc luding the convolu tion of the composite 
spectrum by extinction dPilleri et al.ll2010l) . Figure[3]displays the 
AIB flux, Iaib, obtained by summing the fluxes of the PAH , 
PAH + , and VSG components that were derived from the fit. It 
shows that the AIB intensity correlates well with the [C n] line 
intensity, strongly suggesting that both emissions arise from the 
same regions. The fit of the mid-IR emission provides two inde- 
pendent tracers of the total gas column density N(H) along the 
line of sight as explained below. 

(i) Owing to the excitation mechanism, the AIB intensity can be 
written as \aib k G»x N% ib where N^ /B is the column density of 
carbon in the AIB carriers and Go the UV flux in Habing units 
dHabingj [l968). Assuming that N^ IB /N(H) stays constant at the 
PDR surf ace, Iaib can theref ore be used as a tracer of N(H) if Go 
is known dPilleri et al.ll2010l) . 

(ii) If the column density N(H) is high enough, the effect of ex- 
tinction by silicates can be seen on the AIB spectrum. In the mid- 
IR fit, the extinction is derived from a simple correction term, 
assuming that the emitting and absorbing materials are fully 
mixed: - where tx = N(H) C ext (A) is the optical depth in the 
line of sig ht. The extinction cross-sectio n per nucleon C ext (X) is 
taken from lWeingartner & Draine] d2001l) for R v = 5.5 and N(H) 
is a free parameter of the fit. 

Method (ii) is precise for column densities higher than N(H) ~ 
10 22 cm" 2 . Method (i) can probe lower column densities but has 
two limitations. The AIB emission needs to be corrected for the 
variation in the UV field Go to retrieve the value of N(H). This 
was done assuming that Go scales as the inverse squared dis- 
tance to the illuminati ng star HD200775 and a value of Go =2600 
at 42" from this star dPilleri et al.ll2.Q10t) . We used the projected 
distance as an estimate of the true distance. This introduces an 
error that can be especially strong at positions close to the star 
in the plane of the sky. Figure 3 shows that the AIB emission 
stays almost constant at d < 16", therefore we used this value 
as the minimum effective distance of the NW PDR to the star. 
Method (i) also needs to be calibrated since the local emissivity 
of the AIB carriers is not known precisely. Our approach was 
therefore to derive a calibration factor using the values obtained 
by method (ii) around position 42"on the star-NW cut. The same 
calibration factor was used for all positions along the two cuts. 
Figure|4] shows that the column densities that were derived on 
the two cuts correlate quite well with the [C n] line intensity. 



3. Modelling C + emission 

The critical density for the [Cn] 158 /im line is n e ,„=2500cm~ 3 
for collisions with H and therefore the line emissivity depends 
mainly on temperature for n>n fn/ . We selected a few positions 
on the HIFI S-NW [Cn] cut, three points on the NW PDR and 
two on the South PDR (cf. Table|2]i. The values of Go were deter- 
mined as explained in Sect. 12. 31 and we assumed a constant av- 
erage density with two different values: = 2 x 10 4 cm -3 that 
is characteristic of the molecular cloud dGerin et al ] [19981) and 
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Table 2. Summary of the PDR modelling of the [C n] emission for 5 points along the HIFI S-NW cut (Aa — -47") 



Pos. (AS) 


pro) 


UV field 


T 


PAH+/PAH" 


Ionization 


N(H) 


[C n] local emissivity 




[C n] flux 






(a) 


(b) 


(c) 


parameter (y) (b) 


(d) 


(b) 


HIFI Model (b) 


(") 


(") 


(Go) 


(K) 




(10 3 Go K 1/2 cm 3 ) 


(10 21 cm" 2 ) 


(lO^'WnT 3 ) 


(io- 


7 WnT 2 sr 1 ) 


NW3 (16) 


50 


1873 


337/333 


1.9 


11.0/29.5 


10.5 


4.4/1.4 


8.0 


11.0/8.0 


NW2 (12) 


48 


1975 


342 / 333 


2.7 


11.7/31.1 


8.1 


4.4/1.4 


6.9 


9.5/7.1 


NW1 (-3) 


47 


2100 


348/333 


9.6 


12.6/33.0 


3.9 


4.4/1.4 


3.3 


5.6/4.6 


SI (-63) 


79 


747 


248 / 320 


1.7 


3.8/11.8 


8.7 


4.2/1.4 


3.1 


8.3/7.2 


S2 (-73) 


87 


607 


230/312 


0.96 


3.0/9.5 


14.0 


4.1 / 1.4 


3.7 


9.6/8.4 



(a) Calculated using a projected distance and G =2600 at 42" from the star. 

(b) From the PDR model using n H = 2 x 10 4 / 7 x 10 3 cm 4 , respectively. 

(c) Given as the ratio of the mid-IR intensities shown in Fig.[5] 

(d) Derived from the analysis of the mid-IR emission spectra. 



n H - 7 x 10 3 cm 3 that was derived bv lRapacioli et al. in 
their study of PAH-related species. 

We used the ID Meudon PDR code dLe Petit et alJl2006l) to 
compute the gas temperature T at the cloud surface for all the 
selected positions (cf. Table[2]i. The values of T are used to cal- 
culate the C + level populations. Line intensities are then derived 
by integrating along the line-of-sight (perpendicular to model re- 
sults) and by assuming uniform excitation conditions. The thick- 
ness of the observed regions leads to an optical depth r ~1, 
which implies that transfer effects must be taken into account. 
If we assume constant excitation conditions and gas properties 
along the line of sight, then r and the line intensity can be com- 
puted. 

The agreement between calculated and observed flux values 
is very good when using iih — 7 x 10 3 ctrr 3 . In the NW PDR, the 
ratio is 1.0 for NW3 (16) and NW2 (12), and 1.4 for NW1 (-3). 
For the S PDR, a value of 2.3 is derived for the two positions, 
suggesting that systematic effects are causing the deviation be- 
tween observed and calculated values of the [C n] flux. There are 
several parameters that are not precise in our model but looking 
at Table|2] it seems the local [C n] emissivity is mainly affected 
by the local density and not by the value of Go. For N(H), we 
assumed the same regions emit in PAHs and [C n] , in agreement 
with the profiles shown in Fig. [3] There is also an error for N(H) 
due to our method (cf. Sect. 12. 3t . but this error is expected to be 
the same for both PDRs. Dividing N(H) by a factor of two leads 
to lower values of the ratio of the calculated over the observed 
[CII] flux: 0.7-0.8 for the NW PDR and 1.6-1.7 for the S PDR. 

One step further in the model would consist in studying 
the effect of the grain charge on the photoelectric efficiency 
dBakes & Tielensl il994l) . The relative abundances of PAH + , 
PAH , and evaporating VSGs vary significantly over the neb- 
ula (Fig. [5]). Regions in the cavity appear mainly populated by 
PAH + (cf. NW1 (-3) in Tabled. Since the ionization potential 
of PAH + is much higher than that of PAH (~10eV compared 
to ~6 eV: iMalloci et al.ll2007l) . PAH + should contribute less to 
the photoelectric heating than PAH , leading to a decrease in 
the heating rate, hence in the gas cooling. In its current version, 
the PDR code uses classical grains with an MRN distribution 
dMathis et alj|1977[) and abs orption and scattering cross-sections 
from lLaor & Draind d!993l) . We have used grains of sizes from 
15 A to 3000 A with a dust-to-gas mass ratio of 1%. As a result, 
the ionization parameter y that quantifies the grain charge (cf. 
Table|2]i does not reflect well the variations of the PAH charge 
observed in Fig. [5] An upgraded version, in which the PD R code 
is coupled to the code DUSTEM dCompiegne et al.ll2010l) . is un- 
der development (Gonzalez et al., to be submitted) and will allow 
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Fig. 5. S-NW cut in [Cn] emission measured with HIFI at a 
beam size of 1 1 "(dot-dashed line) and distribution of the emis- 
sion from the different very small dust populations, PAH + (blue), 
PAH (green) and evaporating VSGs (red); filled diamonds are 
Spitzer data (1.8"pix~') and open diamonds are ISOCAM data 
(6"pix- 1 ). 



for including PAHs. NGC 7023 is clearly a template region that 
could be used for these studies. 



4. Conclusion 

By using HIFI and complementary mid-IR data, we have shown 
that the [C n] cooling line and the AIB emission arise from the 
same regions, in the transition zone between atomic and molec- 
ular gas. The prototype PDR NGC 7023 was found to be a good 
object for comparison with PDR models. Further progress on the 
energetics of this region awaits for the coming [O i] data from the 
PACS instrument and a PDR model that treats the photophysics 
of PAHs consistently. 
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